273 research outputs found
IFCC recommendation on sampling, transport and storage for the determination of the concentration of ionized calcium in whole blood, plasma and serum
The substance concentration of ionized calcium (c
Ca2+) in blood, plasma or serum preanalytically may be affected by pH changes of the sample, calcium binding by heparin, and dilution by the anticoagulant solution
Coupled coarse graining and Markov Chain Monte Carlo for lattice systems
We propose an efficient Markov Chain Monte Carlo method for sampling
equilibrium distributions for stochastic lattice models, capable of handling
correctly long and short-range particle interactions. The proposed method is a
Metropolis-type algorithm with the proposal probability transition matrix based
on the coarse-grained approximating measures introduced in a series of works of
M. Katsoulakis, A. Majda, D. Vlachos and P. Plechac, L. Rey-Bellet and
D.Tsagkarogiannis,. We prove that the proposed algorithm reduces the
computational cost due to energy differences and has comparable mixing
properties with the classical microscopic Metropolis algorithm, controlled by
the level of coarsening and reconstruction procedure. The properties and
effectiveness of the algorithm are demonstrated with an exactly solvable
example of a one dimensional Ising-type model, comparing efficiency of the
single spin-flip Metropolis dynamics and the proposed coupled Metropolis
algorithm.Comment: 20 pages, 4 figure
Conformational and Structural Relaxations of Poly(ethylene oxide) and Poly(propylene oxide) Melts: Molecular Dynamics Study of Spatial Heterogeneity, Cooperativity, and Correlated Forward-Backward Motion
Performing molecular dynamics simulations for all-atom models, we
characterize the conformational and structural relaxations of poly(ethylene
oxide) and poly(propylene oxide) melts. The temperature dependence of these
relaxation processes deviates from an Arrhenius law for both polymers. We
demonstrate that mode-coupling theory captures some aspects of the glassy
slowdown, but it does not enable a complete explanation of the dynamical
behavior. When the temperature is decreased, spatially heterogeneous and
cooperative translational dynamics are found to become more important for the
structural relaxation. Moreover, the transitions between the conformational
states cease to obey Poisson statistics. In particular, we show that, at
sufficiently low temperatures, correlated forward-backward motion is an
important aspect of the conformational relaxation, leading to strongly
nonexponential distributions for the waiting times of the dihedrals in the
various conformational statesComment: 13 pages, 13 figure
Molecular Dynamics Simulations
A tutorial introduction to the technique of Molecular Dynamics (MD) is given,
and some characteristic examples of applications are described. The purpose and
scope of these simulations and the relation to other simulation methods is
discussed, and the basic MD algorithms are described. The sampling of intensive
variables (temperature T, pressure p) in runs carried out in the microcanonical
(NVE) ensemble (N= particle number, V = volume, E = energy) is discussed, as
well as the realization of other ensembles (e.g. the NVT ensemble). For a
typical application example, molten SiO2, the estimation of various transport
coefficients (self-diffusion constants, viscosity, thermal conductivity) is
discussed. As an example of Non-Equilibrium Molecular Dynamics (NEMD), a study
of a glass-forming polymer melt under shear is mentioned.Comment: 38 pages, 11 figures, to appear in J. Phys.: Condens. Matte
Simulation of dimensionality effects in thermal transport
The discovery of nanostructures and the development of growth and fabrication
techniques of one- and two-dimensional materials provide the possibility to
probe experimentally heat transport in low-dimensional systems. Nevertheless
measuring the thermal conductivity of these systems is extremely challenging
and subject to large uncertainties, thus hindering the chance for a direct
comparison between experiments and statistical physics models. Atomistic
simulations of realistic nanostructures provide the ideal bridge between
abstract models and experiments. After briefly introducing the state of the art
of heat transport measurement in nanostructures, and numerical techniques to
simulate realistic systems at atomistic level, we review the contribution of
lattice dynamics and molecular dynamics simulation to understanding nanoscale
thermal transport in systems with reduced dimensionality. We focus on the
effect of dimensionality in determining the phononic properties of carbon and
semiconducting nanostructures, specifically considering the cases of carbon
nanotubes, graphene and of silicon nanowires and ultra-thin membranes,
underlying analogies and differences with abstract lattice models.Comment: 30 pages, 21 figures. Review paper, to appear in the Springer Lecture
Notes in Physics volume "Thermal transport in low dimensions: from
statistical physics to nanoscale heat transfer" (S. Lepri ed.
Combining Optimal Control Theory and Molecular Dynamics for Protein Folding
A new method to develop low-energy folding routes for proteins is presented. The novel aspect of the proposed approach is the synergistic use of optimal control theory with Molecular Dynamics (MD). In the first step of the method, optimal control theory is employed to compute the force field and the optimal folding trajectory for the atoms of a Coarse-Grained (CG) protein model. The solution of this CG optimization provides an harmonic approximation of the true potential energy surface around the native state. In the next step CG optimization guides the MD simulation by specifying the optimal target positions for the atoms. In turn, MD simulation provides an all-atom conformation whose positions match closely the reference target positions determined by CG optimization. This is accomplished by Targeted Molecular Dynamics (TMD) which uses a bias potential or harmonic restraint in addition to the usual MD potential. Folding is a dynamical process and as such residues make different contacts during the course of folding. Therefore CG optimization has to be reinitialized and repeated over time to accomodate these important changes. At each sampled folding time, the active contacts among the residues are recalculated based on the all-atom conformation obtained from MD. Using the new set of contacts, the CG potential is updated and the CG optimal trajectory for the atoms is recomputed. This is followed by MD. Implementation of this repetitive CG optimization - MD simulation cycle generates the folding trajectory. Simulations on a model protein Villin demonstrate the utility of the method. Since the method is founded on the general tools of optimal control theory and MD without any restrictions, it is widely applicable to other systems. It can be easily implemented with available MD software packages
Systematic Computational and Experimental Investigation of Lithium-Ion Transport Mechanisms in Polyester-Based Polymer Electrolytes
Understanding the mechanisms of lithium-ion transport in polymers is crucial for the design of polymer electrolytes. We combine modular synthesis, electrochemical characterization, and molecular simulation to investigate lithium-ion transport in a new family of polyester-based polymers and in poly(ethylene oxide) (PEO). Theoretical predictions of glass-transition temperatures and ionic conductivities in the polymers agree well with experimental measurements. Interestingly, both the experiments and simulations indicate that the ionic conductivity of PEO, relative to the polyesters, is far higher than would be expected from its relative glass-transition temperature. The simulations reveal that diffusion of the lithium cations in the polyesters proceeds via a different mechanism than in PEO, and analysis of the distribution of available cation solvation sites in the various polymers provides a novel and intuitive way to explain the experimentally observed ionic conductivities. This work provides a platform for the evaluation and prediction of ionic conductivities in polymer electrolyte materials
- …